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Abstract 



The chiral phase transition of the finite temperature QCD is discussed in the lad- 

P-i! 

der approximated Exact (Wilson's) Renormalization Group (ERG). The gluon thermal 
mass is incorporated by the Hard Thermal Loop (HTL) approximation. We calculated 
a critical temperature T c , a critical exponent u, and order parameters such as chiral 
condensate < ipip > and pion decay constants with this approximation. The depen- 
dences of the pion decay constants on an unphysical parameter are fairly improved by 
including contributions from momentum dependent vertices. The critical temperature 
was found as: T c w 124 MeV for 3-flavors. It was shown that the critical tempera- 
ture considerably depends on the approximation scheme for the gluon thermal mass 
term. We also calculated the temperature dependence of the chiral condensate, 7r decay 
constants, and the pion velocity. 
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§1. Introduction 



Explorations of the chiral nature of QCD at high temperature are experimental as well as 
theoretical interests. They are relevant to the early universe and to the heavy ion collisions 
at RHIC and those planned at LHC. It is believed that the approximate chiral symmetry 
observed at bare Lagrangian dynamically breaks down at zero temperature due to the chiral 
condensate contained in the QCD vacuum. At high temperature, the chiral condensate is 
expected to melt, and the chiral symmetry is restored at a certain critical temperature. Since 
the perturbative methods can not be applied to the non-trivial QCD vacuum, so we need 
non-perturbative methods, such as the lattice Monte-Carlo simulation, the Schwinger- Dyson 
equation, and the Exact Renormalization Group, for investigating the vacuum structure of 
QCD at high temperature. 

For the zero temperature case, the Schwinger-Dyson equations (SDEs) has been stud- 
ied 1 )" 4 ) and has been applied to the strong coupled QED and QCD. 5 )" 7 ) Commonly, the 
so-called (improved) ladder approximation 5 ) has been applied to solve the SDEs. Although 
these approximation schemes seem to be fairly crude, the SDEs offer surprisingly good re- 
sults on the dynamical chiral symmetry breaking in QCD. For the finite temperature gauge 
theories, the SDEs were also applied. 8 )" 14 ) At finite temperature, the gauge interaction is 
corrected by thermal effects. Thermal effects, such as the Debye screening and the Landau 
damping, appear in the gauge boson propagator as their thermal mass. So it is no longer 
to say that how to approximate the thermal mass is one of the main subjects for the im- 
provement of finite temperature analyses. The hot QED was discussed in the instantaneous- 
exchange (IE) approximation for the gauge boson propagator including the screening effect in 
the static limit, 8 ) or the IE approximation for the longitudinal (electric) mode and the Hard 
Thermal Loop Approximation for transverse (magnetic) modes. 9 ) The chiral phase transi- 
tion of the hot QCD was studied by neglecting the gluon thermal/Debye mass term. 10 )" 12 ) 
The Phase structure of hot and/or dense QCD was explored by introducing the Debye mass 
in static limit of the gluon propagator. 13 ) However unfortunately the improvement of the 
above approximations seems to be difficult due to the complexity of the integral equation. 

The Exact Renormalization Group (ERG) has been applied to investigate the chiral 
nature of the strong coupling gauge theory at zero temperature. 15 ) With ERG method, 
we can study the chiral critical behavior in very simple manner, and if we approximate 
the ERG flow equations to the ladder-like corrections then it leads the identical results 
with those obtained by solving the ladder SDEs. It should be noted that we can improve 
the approximation scheme systematically, e.g. incorporating the "crossed-ladder" diagrams. 
Especially, the large gauge dependence of the solutions obtained by the ladder SDEs 16 ) has 
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been improved in the case of the chiral critical behavior in QED. 17 ^ 

For the finite temperature, the chiral nature has been investigated by the "Quark-Meson 
model" in ERG method. 18 ) The Quark-Meson model is defined at 'so-called' compositeness 
scale kq, « 600 MeV from the zero temperature QCD and therefore does not incorporate 
the thermal screening of the gluon interaction. It is expected that the mesonic bound states 
forms at the short scale region: k > kq> and are not affected by the long distance: k < kq> 
behavior of the gluon propagator. It is one of the issue whether the thermal screening effect 
to the chiral behavior can be neglected or not. 

In this paper, we investigate the chiral phase transition of the hot QCD in the ladder 
approximation with including the gluon thermal mass in the Hard Thermal Loop (HTL) 
approximation. To introduce auxiliary fields, we consider a intermediate momentum scale 
similar to k$ in the Quark-Meson model. The pion decay constants depend on the interme- 
diate scale, which is an unphysical parameter. By including contributions from momentum 
dependent vertices, the decay constants become considerably stable for the intermediate 
scale. The formation of the chiral condensate is remarkably affected by the thermal screen- 
ing. Especially, the critical temperature is remarkably lowered in comparison with that 
given by neglecting the thermal mass. We also calculated the temperature dependence of 
the order parameters, i.e. the chiral condensate and the pion decay constants for several 
approximation schemes for the thermal mass. 

We estimated the temperature dependence of the pion velocity too, from the ratio of 
the spatial component of the decay constant and the temporal component of that, which 
become different at finite temperature. 19 - ) By the ladder approximated SDE, the temper- 
ature dependence of the 'space-time averaged' component of the pion decay constant was 
evaluated, 11 ^ however, that of the pion velocity was not. It is expected that pions travel 
at less than the speed of light in a thermal bath. In this paper, it is shown that pions are 
slowed down slightly near T c . 

This paper is organized as follows. We explain our approximation scheme and derive 
the ERG flow equations in £0 in SI we show the numerical results. 21 is devoted to the 
summary and discussions. 

§2. The ERG equation 

In this section we explain the approximation scheme first, and next we derive the ERG 
equation. Let us start form iVy-flavor QCD at finite temperature T. The bare Lagrangian 
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density C is given by 

C = \tvF^ + qiJpq + i- (cyi^) 2 + [d,5 AB + gf ABC AC] C B , (2-1) 

where q denotes the Nf x N c components quarks. N c is the number of color. The field 
strength F^ u = F^ V T A and the covariant derivative are given by, 

D, = d, + igA*T A , (2-2) 
K = d ^ ~ d »K + 9fAB C A°A*. (2-3) 

The third term is gauge fixing one, and the fourth is from FP ghosts. We choose the Landau 
gauge a = 0. Here C is invariant under the chiral flavor symmetry Ur{Nj) x U^Nf)*' . The 
axial current is defined as J5 = qi/y^q, where a denotes iso-spin indices and A a are the 
generators of the chiral U(N/) (trA a A fe = \5 ab ). 

2.1. Gauge interaction 

There are some problems in the application of the ERG to the gauge theory because 
of the infrared cutoff term which explicitly breaks gauge invariance. 20 -'' 21 -' There are two 
ways to treat the gauge theory. One is generalizing the ERG keeping the gauge symmetry. 21 ) 
This is of course an ideal method, however their formulations are not accompanied with the 
non-perturbative approximation method. The other is introducing the gauge non-invariant 
counter terms to compensate the gauge invariance. Hence, the bare action does not satisfy 
the Slavnov- Taylor Identity, but Modified one. 20 ) However, since we would like to introduce 
the gluon self energy in the HTL approximation, we must calculate the momentum depen- 
dent gluon self energy and therefore relevant terms of the Modified Slavnov- Taylor Identity 
(MSTI) become complicated integral equations. 

In this paper, we do not adopt the either methods explained above. We integrate out 
gluons at the beginning and solve ERG flow equations for the reduced fermionic theory. This 
procedure generally causes many multi-fermi interactions. By the diagrammatical consid- 
eration, one can easily realize that all ladder type diagrams are generated by only the four 
fermi operator induced by the one gluon exchange. Other multi fermi operators generate the 
non-ladder diagrams. So as shown in Fig. Q we replace the gauge sector with the non-local 
four fermi interaction induced by the one gluon exchange i.e. 

Gi = I I I g(p,k)g(q,-q-p-k)q(k) 1 ,TMp)K B (P + k )^^ ( 2 ' 4 ) 

J q J p J k 

*' The axial chiral symmetry U(1)a is broken due to the anomaly, and the QCD phase transition is 
affected by the instanton effect. However, we do not consider the instanton effect in this paper. 
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Fig. 1. 1-gluon exchange as 4-fermi interaction in the ladder approximation. 



where g(p, k) and are the running gauge coupling constant and the gluon propagator 
respectively. Note that, since other corrections from the gluon fluctuations do not contribute 
to the ladder approximation, our treatment does not mean the additional approximation. 
Now, what we should investigate is reduced to the fermionic theory interacting through the 
non-local four fermi interaction ()2-4j) and the ladder approximate solutions can be found in 
the leading order of 1/N C expansion with the coupling constant g: g 2 = g 2 /N c . 

The flow of the momentum dependent four fermi interaction is studied by Meggiolaro 
et al. 22 -* They do not integrate out gluons at the beginning. They introduce the IR cutoff 
term for gluons and the ERG flow starts from the QCD classical action. In their work, to 
compensate the gauge invariance, the gluon mass is given so as to satisfy the MSTI at the 
leading perturbative order. We do not examine the flow of such interactions themselves as 
they did, but the flow of the effective potential V. The contribution from the momentum 
dependent four fermi interaction appears in the ERG equation for V and affects the chiral 
transition, in which we are interested. 

The argument of the gauge coupling constant should be taken as g(p, k) = g((p + k) 2 ) to 
satisfy the chiral Ward identities. 23 -* The infrared singularity of the gauge coupling constant 
is regularized as : 



9 2 (A) 



61n(AV^ CD ) if A > Ax 



1 



+ \^A,iA^tA^A-iJ if A X >A>A 2 (2-5) 



61ii(^iMqcd) 1 61n(^ 2 /A 1 )(ln(^ M 2 jCD ))2 

1 HM/Ai) if A < A 

61n(^M* CD ) b(HAi/A- QCD )r ~ /l2 



where ln(yVvl QCD ) = 0.25, ln(A 2 /A QCD ) = -1.0 and b = (11N C - 2N f )/48** 
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At high temperature gluons become massive by the thermal fluctuation. We incorporate 
the Debye-screening effect to the gluon propagator D£„ in the HTL approximation. Of 
course, the HTL approximation is based on the perturbative QCD and it is valid in the high 
energy (short distance) region. However, it is expectable that the short distance behavior 
is important to the QCD chiral physics, not the long distance one. Also, it is known that 
the order parameters, such as the chiral condensate, hardly depend on the difference of the 
treatments of the running coupling constant in the low energy region. 7 ) So applying the HTL 
approximation for the gluon propagator is a proper procedure for our present purpose , which 
is the first step of the improvement of the approximation for introducing thermal effects. Of 
course it should be checked that how order parameters and the critical temperature depend 
on the IR cutoff scheme for g(A) at finite T. We show the Ai dependence of our numerical 
results later, in 

The gluon self-energy is a gauge dependent quantity, contrary to the photon self-energy 
in QED. However the HTL approximated one has no dependence on the choice of gauge. At 
finite temperature, the self-energy matrix takes the following form: 

n£={FI% + GI%)6 AB . (2-6) 

The projectors P^ v and P^ u are given by 

P T — A — Pj^j pL _ X _ PfJ-Pv _ pT Icy y\ 

r fjiv u fMV - 2 > r iiv u in> ^2 \a>' \*'*) 

We use the following notation, p^ = (0,p) and = diag(0, 1, 1, 1). One can easily verify 
that = 2, P^ = 1. In the HTL approximation, F and G take the form: 24 ) 



F = p 2 + ^n 00 (p ,p) =p 2 + ml Th (rj 2 - 1)<P( V ), (2- 



G = p 2 + ~n u (p ,p) - jii7 00 (po,p) = P 2 + \rn 2 HTL [l - ( V 2 - 1)$( V )] , (2-9) 
where 77 = ipo/\p\. At temperature T, the thermal mass m HXL and $(rj) are given by 

< TL = \g 2 T 2 (n c + ±N f ) , (2-10) 



3 V 2 
^) = ^ln 



2 \7] - 1 

fUan-M.!. (241) 
\P\ Po 
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2.2. The auxiliary fields 

It is convenient to introduce the auxiliary fields for the quark- ant iquark pairs to make cal- 
culation efficient. 15 ) We introduce the auxiliary fields for the composite operators <Lt AUX <Z J d AUX 
and iqA Aux j5qA A uxj where qA AVX is the low energy mode of the quarks defined by 

T E / J^(pK^aux -P). (2-12) 

The q integral reduces to the Matsubara sum, i.e. q — > 2{n + 1)ttT, n = 0, ±1,±2, 
The intermediate momentum scale A AVX plays similar role to the compositeness scale /c$ in 
Ref.18). However we do not impose the compositeness conditions 25 ^ at this scale. 

The auxiliary fields a and 7r are introduced by inserting the following trivial Gaussian 
integrals: 

l = N ~ X j ^expj-^M 2 J^dr j d 3 x(a-yM- 2 q AAUX q AAVX ) 2 y (2-13) 
1 = * f ~ 1 J Dnexp^M 2 J^dr J d 3 x (tt - yM~ 2 iq AAVxlb q AAVX ) 2 ^ , (2-14) 

where J\f is a irrelevant constant and (3 = Since the compositeness conditions are 

not imposed, M 2 and y are the arbitrary constants. Of course, the physical quantities should 
be independent of M 2 , y and A AUX . We introduced the composite operators whose tensor 
structures are 1 and 75, and neglected the others, like the Pagels-Stokar approximation. 26 > 
Using G\ of (|2-4j) and the contribution from (|2T3J) and (j2-14j) . we can write our Lagrangian 
density C A as 

£ A = qi$q + V a + G 1 , (2-15) 

where 

V a = \M- 2 y 2 {{q AAVX q AAVX f + (iqA AVX 7 5 qA AVX ) 2 ) 

-y{<7QA AVX qA AVX + -xiq AAVX ^q AAVX ) + ^M 2 (a 2 + n 2 ). (2-16) 

The ERG evolution is divided into two regions, A > A AVX and A < A AVX . In the high 
energy region A > A AVX , C A describes a pure fermionic system interacting through G\. At 
A = A AX jx, additional Yukawa interactions and four fermi interactions are switched on. In 
A < A AUX , our model becomes the quark-meson system with G\. We separate quark fields 
into the high energy mode and the low energy mode at the scale A AVX , and the auxiliary 



*' In the following, we choose units such that Boltzmann constant k = 1. 
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fields for mesons are introduced as what is coupled with the low energy mode, ?a A ux- We 
show the A Avyi dependence of the order parameters in 

The intermediate scale A AVX plays almost similar role for k<$> in the Quark-Meson model, 18 ) 
however A AUX is not an UV cutoff for some effective theory. We solve the ERG flow equations 
from the sufficiently large cutoff A ^> yl AUX for each temperature T. In both regions: 
p > A AUX and p < A AUX , the gluonic contributions are explicitly evaluated through one gluon 
exchange term G±, since we are interested in the influence of the gluon thermal screening 
effect for the chiral phase transition. 

2.3. Flow equations in the ladder approximation 

Now, let us derive the ERG equation. The flow equation describes the change of the 
effective average action r A under an infinitesimal shift of the infrared momentum cutoff A. 
We start with the construction of r A . The generating functional of the connected Green 
function is 

W A [J] = \nj DqDqexp{-S cut - S barc + J dr J d 3 x(fjq - q V )}, (2-17) 

where Share is the bare action, i.e. Share = / drd 3 x\CA, and {fj, 77} are external sources. The 
IR cutoff term S CVLt is introduced for the (anti-) quark fields, 

Scut[q,q} = J^dr J rf 3 x qA-\-id, A)q. (2-18) 

Here A -1 is a chiral invariant cutoff operator and has the property, 

* i, .s f for p^>A , 
A-\p,A)^\ P (2-19) 

I 00 for p A. 

We do not introduce the IR cutoff term for boson fields in this paper since we treat the 
ladder part only. In the ladder approximation, the bosonic degrees of freedom, a and tt are 
not integrated in (J2-17|) . So we can forget the infrared problem for mesons at T < T c . The 
infrared problem originated from the absence of the IR cutoff for gluons can be regularized 
by introducing a sufficiently small IR cutoff Ai& to the gluon propagator in G\. However, 
our results are not affected by /Iir, since our ERG evolutions of the effective potential as 
well as the order parameters are frozen fast when A gets smaller than the typical scale for 
the chiral physics. 

The cutoff effective action r& is defined by Legendre transformation 

r A [0] = -W A [J] +J P drJ rf 3 xJ • $ + S cut [<P], (2-20) 
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where <P is given by SWa[J]/SJ and we use the shorthand notation of the external sources: 
J = {f],r] T }- Note that, since >S cu t preserves the chiral symmetry, the effective action also 
respects it. 

Differentiating both sides of Eq. (j2-17j) and using (j2-20|) . we get the ERG equation for the 
effective average action, 27 ) 



A±r 4 *] = -istr 



,(2) 



(2-21) 



where Str is super-trace which involves momentum (or coordinate) integration, Matsubara 
summation, spinor summation and color summation. r\ is a second (functional) derivative 
with respect to the fields <P T = (q T , q), i.e. 



r 



(2) 



■'■y 



5$,. 



(2-22) 



We employ the sharp cutoff scheme in this paper since in the smooth cutoff case we 
must perform the momentum integration with respect to the spatial momenta numerically 
for each Matsubara mode. We want to avoid such a tiresome and unpractical work. The 
sharp cutoff limit of Eq. (|2-2I|) can be written as: 28 ) 



~dA 



A 



Str 



8J\q\-_A) 

7(q) 



(2-23) 



where is a interaction part of the second derivative of and 7(q) is an inverse propa- 
gator: 

' o i 

e o 



7(q) 



(2-24) 



The infrared cutoff propagator C(q) is given by, 

9{\q\-A) 



C(q) 



7(q) 



(2-25) 



It is known that, in the sharp cutoff case the vertices of the effective action i -1 ^ should be 
expanded in powers of \p\ = a/p 2 instead of the derivatives c^. 28 -* Heaviside 9 function in 
ERG equations is also expanded in terms of |p| as 



9(\p -q\-A) = 9(p 2 + 2p-q) = 9(p-q) + J2 ^ n) {p ■ q) 



(2-26) 



where p is a unit vector parallel to p. Integrating r.h.s. of the ERG equation ()2-23J) with 
respect to the internal momenta q, one can expand it in terms of the momentum scale \p\. 
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We approximate the effective action as: 

r A = J dr J d 3 x {qi$q + V(q, q, a, 7T ) + V a + G 1 + K(ir, iqj 5 q)} , (2-27) 

where K(ir, iq^q) is the higher order term of momentum scale expansion. The explicit form 
of K is shown later, in ()2-38|) . We neglect the kinetic term of a, which is not concerned 
with the analysis here. The potential part V(q, q, a, it) is a function not only of S = qq and 
P = «975<7 but also of ?yi AUX ?4 AUX and iqA AUX 75qA AUX , therefore vertices of V have momentum 
dependence like 6*(yl AUX — p). 

The complexity arising from 8(A AVX — p) can be avoided by dividing the ERG evolution 
into two regions, i.e. the high energy region A > A AVX and the low energy region A < A AVX . In 
the high energy region, <?/i AUX <?/i AUX and «<Ld AUX 759A AU x d° not contribute to r.h.s. of Eq. (J2-23|) 
since the supports of q~A AUX (p) and q AAVX (p) are restricted to p < A AVX < A. Thus we can 
forget V a in this region. Next, in the low energy region, we need not distinguish qA AVX (p) 
from q(p) in the vertices of the effective action with external momentums Pi < A AX]X and of 
course in the ERG equation (|2-23j) . Hence we can replace qA AVX (p) with q(p) by choosing 
A AVX larger than typical scales of the chiral dynamics. Now our prescription is that in the 
high energy region we solve ERG equation by neglecting V a term and at A = A AVX we shift 
the potential V(q, q, a, n) — > V(q,q,a,7r) + V a (qj\ AVX = <?, 1a avx = <?)• Hereafter in both 
regions, we write the arguments of the potential V as S, P, a and 7r i.e. V(S, P, a, n) 

We can get the ERG equation for V by substituting the zero modes of the composite 
operators S = J drd 3 xS(x) and P = J drd 3 xP(x)*\ The ERG equation for V can be 
read, 28 ) 

A lA V = "i Str [ ln ( 7(q) + V{2) + G ^)] ' (2,28) 

where and are the matrices of the second derivative of V and of the four fermi 
interaction induced by one gluon exchange ()2-4|) with respect to the (anti-)quark fields re- 
spectively. 

The derivatives of V with respect to (anti-) quarks can be written as 

' ] - V ^n = Wj4t;V(S, P) + i^jl^V^ P) + ~-, (2-29) 



dqi, a dq J > b b J dS K ' ' b J '°dP 



*' For the fermionic theory, we can not insert the zero modes of the fermion fields. Otherwise, the higher 
vertices which is necessary to produce the bulk structure of the general Green function of the fermions vanish 
due to the Grassmann nature of fermions. 
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where we omit the terms generating the non-ladder diagrams in ERG. Thus the matrix T/ (2) 
is 

V W = ( WAVs + vysVp)} 

\ Sffl (V s + tilV P ) J ' 

where subscripts 'S" and 'P' denote the derivative with respect to S = qq and P = iq^^q 
respectively, e.g. V$ — dV/dS. The second functional derivative of the gauge induced four 

(2) 

fermi interaction G\ is calculated as: 

The argument of the gauge coupling constant g is an I.R. cutoff A, i.e. leading order in the 
momentum scale expansion. The ladder contributions to a channel and tt channel are given 
as follows, 

-T A lu q J ® q jlfl T A = --J—S^S'jS^qq + ij 5 qij 5 q) l c + • • • , (2-33) 



where we use the relation of the generator T . 

(T A )i(T A ) l k = - ±-5{5l, (2-34) 

and l c is a color unit matrix. Consequently, we find the ladder contribution from the gauge 
interaction G\ as: 

r (2) _ _ 9 2 (A) ( 5&(S + * 7s P) \ 

1 " 8iV?JV> ^ )aa V -W(5- + i T5 T P) )■ { } 

The trace of the gluon propagator tr_D _1 = {D~ X Y/ A reduces to a factor 

F G A 2 + m|(p /|p|) /l 2 + m 2 (p /|p|) 1 J 

for iV c = 3, where and my are longitudinal component and transverse component of 
thermal mass respectively. The argument Po/[p| of thermal mass is written in terms of 
A and u n = (2n + l)ixT : (^/ A 2 /uj 2 — 1)~ 1//2 in the ERG flow equation. The bosonic 
Matsubara frequency is u n = InixT . However, in our analysis, the momentum which flows 
into G\ in (|2-4jl is expected to (2n + i)irT because of the derivative expansion. We expand 
Gi((2n + 1)tvT + Ui) in terms of ui, which is odd, and include only the leading part: I = 0. 
Inserting Eqs. flZZS) , d2~TTT|l and (|2~35|) to Ea. <|72B)l . the final expression of the potential 
part of the ERG equation reduces to 

= ^/^E'v^^I 1 * + Vs + m , (2-37) 
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where we write V = V — g 2 (S 2 + P^trD^ 1 /32N 2 N f and J^ n ' is a summation over Mastubara 
modes n with the condition A 2 > u 2 . Thus the RG evolution terminates at A — ttT. 

The pion wave function renormalization factors of the spatial derivative and the temporal 
one are different because the Lorentz invariance is lost. For the derivation of the ERG flow 
equations for the 7r wave function renormalization factors, Z^l and Z„t, we consider the 
higher order contribution of momentum scale expansion in ()2-27j) . 

K(tt, P) = iz w (p)7r p 7r_ p + K(p)7r p P_ p + ^G(p)P p P^ p , (2-38) 

where P p is a pionic composite operator: P p = (qi'j5q) p , and subscripts p and — p denotes 
their momenta. Y v (p) and G(p) are the momentum dependent vertices of Yukawa interaction 
and that of four fermi interaction respectively. Taking functional derivative of K(tt, P) with 
respect to P p and P- p , and substituting to ()2-23|) . the ERG equation for K is given by 



A JA K(P) = | St 4 5 (l^|-^) S ^(q)^(b+g|-^)Si,(p + q)Y]^^(yl AUX -yl). (2-39) 

K denotes the 'interaction parts' of K, where the pion kinetic term is subtracted in ()2-38|) . 
The step function 9(A AUX — A) in r.h.s. of Eq. (j2-39|) is due to the cutoff of the quark-meson 
interaction. It tells that K contributes to the whole ERG flow for A < A Avyi . Sf and Y 
are the massive quark propagator on the vacuum < a >= oo and the matrix part of the 
Yukawa coupling constant: Y n (p) = Y(p)Y respectively. Introducing the coupling constants 
V = m qq + ■ ■ ■ , is written as: 




S/(q)= 7 (q)+ „ , (2-40) 



and Y is given by, 





T 




(2-41) 



To devive the ERG equation of Z n i and Z n x, w e expand (|2-39|) in powers of momentum 
scale to second order. Note that, terms which mix p and p do not contribute to the wave 
function renormalization factors. The momentum dependent vertices are expanded as: 

Z K (p) = Z wL p 2 + Z nT p 2 + 0(p), (2-42) 
Y(p) = (\ Po \/A)y lL + (p/A)y lT + (p /A) 2 y 2L + (p/A) 2 y 2T , (2-43) 
G(p) = (\p Q \/A)g 1L + (p/A)g 1T + (p /A) 2 g 2L + (p/A) 2 g 2T . (2-44) 
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We need to be careful about expanding the inside of square bracket in ()2-39|) in terms of po, 
which takes discontinuous value at finite temperature. Let us expand a function f(po) as: 



/(Po) = /o + /ibo| + / 2 Po- (2-45) 

We estimate the expansion of f(po) using three points of bosonic Matsubara frequency, 
Po = 0, po = 2ttT and po = 47rT, as: 

47TTA = -3/(po = 0) + 4/(po = 2ttT) - f(p = 4ttT), (2-46) 
8n 2 T 2 f 2 = f(p = 0) - 2/(po = 2ttT) + f(p = 4ttT). (2-47) 

The result of the expansion above becomes identical to that of the ordinary momentum scale 
expansion at T — > 0. Finally, we get following ERG equations, 

A ~IA ZlT{L ' T) = m2/2 ( z " T ) + ( 2m ?/2(L,r) + 2/i(l,t)) j o + 2rnynL,T)h(L,T), (2-48) 

A jA yi(L ' T) = m9aIl{ - L > T ) + (y^L,T)9a + mg 1{LiT) )I 0l (2-49) 

A JA V2{L,T) = mg<jI ^ L ^ T ) + (Vi(L,T)g<, + mgnL tT ))h(LX) 

+(rng 2 (L,T) + Vi{l,T)9i{l,T) + y2(L,T)9a)Io, (2-50) 

^TT^i(L,T) = glh(L,T) + 2g*9x(.L,T)Io, (2-51) 



A dA 92(L,T) = 9 ° l2(L ' T) + 2 9^9\{L.T)h{L.T) + (2g a g 2 ( L ,T) + gl( L ,T))h- (2-52) 



The threshold functions Iml,t) include Matsubara sum. The explicit forms of Ik(L,T) are 
given in Appendix El ERG equations for m and g a are derived by the polynomial expansion 
of V in terms of the composite operator. We show the explicit form later. 

§3. Numerical Analysis 

3.1. Truncation and A AUX dependence of the order parameters 

Before starting the analysis of the chiral phase transition at finite temperature, we in- 
vestigate the truncation and A AVX dependence of the order parameters. In general, the 
approximate solutions depend A AUX since the A AVX dependent factors can not be eliminated 
by the normalization of the redundant degrees of freedom a, tt. Hence the ^l A ux dependence 
of the results should be verified. In this paper we show the truncation//l AUX dependence of 
the order parameters such as the pion decay constant /„- and the chiral condensate < ipip > 
at zero temperature. 
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At zero temperature, the ERG flow equation for V reduces to 

A±V=±N,N e la(A> + V§ + V2), (3-1) 

where V = V — 3g 2 (S 2 + P 2 )/2N 2 NfA 2 . This flow equation leads to the ladder-exact results 
for the chiral condensate < ipip > and the dynamical quark mass m e j. 15 ' Here (?/>) ip denotes 
the (anti-) quark field with a definite flavor. 

To estimate the ir decay constant, we must calculate the wave function renormalization 
factor of the meson fields Z n . The ERG equations of Z n and momentum dependent vertices 
are the same form as ()2-48j) - ()2-52j) . At zero temperature, there exists no distinction between 
L and T. The explicit forms of zero temperature coefficients Ik are also given in Appendix 

El 

The 7r decay constant f n is given in terms of Z n and the vacuum expectation value (VEV) 

f w = y/2/3Z^ 2 a . (3-2) 

The coefficient differs from that in 18) because our normalizations of the auxiliary fields are 
different. To estimate the chiral condensate, we introduce the external source mo for the 
composite operator qq. The chiral condensate is given by 

- < >= -J, «£W , (3-3) 

Nf om Q m =o 

where we write the renormalization point A since this formula leads the chiral condensate 
renormalized at U.V. cutoff scale A . The ERG flow starts from A = A , so the solutions 
of the ERG depend on A . In this paper, we set A as: A = A QCD x e 10 ' 5 . Hereafter, we 
consider the case: Nf = 3. 

We first investigate the truncation dependence of the order parameters. The potential V 
is truncated to the (N — l)-th polynomial of S, 

N 

V(S) = Y,v n -iS n -\ (3-4) 

n=l 

where v n are the VEV cxq dependent coefficients. By virtue of the parity invariance, the 

higher terms of P do not contribute to 0(P°) terms. So we can neglect such terms in (|3-4|) 

without losing generality. The ERG flow equation reduces to a coupled ordinary differential 

equation. The truncation (N) dependence of f n is shown in FigO As shown in Fig|21 tt 

decay constant rapidly converges against truncation of the effective potential V(S). We set 

*' We introduced flavor-singlet auxiliary fields by (|2-13|) and (|2-14|) . In our approximation, the decay 
constant can be regarded as that of meson flavor-5't/(3) octet. 
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^Iqcd to 583 MeV, so as to the convergent f n becomes 92.42 MeV, which is the experimental 
value at T = 0. 29 ) Throughout this paper, yl QC D is fixed to this value. In this calculation, 
we choose At = ln(yi AUX //lQCD) = 1-0. 

140 

130 

> 

J 120 
^ 110 
I 100 

« 

§ 90 
o 

> 80 

to 

o 

& 70 

o 60 
a 

50 
40 

3456789 10 
Order of Truncation (N) 

Fig. 2. Truncation dependence of the it decay constant. The dotted line is f w = 92.42 MeV. The 
order of truncation corresponds to N. 

In TableHJ the truncation dependence of the chiral condensate, the quark dynamical mass 
and the it decay constant are described. The chiral condensate is renormalized at 1 GeV. 
Although the chiral condensate and the dynamical quark mass have already been calculated 
in Ref.15), we show the truncation dependence of these quantities since the results may 
depend on the unphysical parameters, e.g. At, for the lower truncation N < 5 where the 
results do not converge yet. 

Table I. Truncation dependence of the chiral condensate, the quark dynamical mass and the ir 
decay constant, if) denotes the quark with a definite flavour. 



Truncation (N) 


4 


5 


6 


7 


8 


9 


10 


| < ^ > | 1/3 (MeV) 


273.2 


258.2 


260.8 


259.6 


260.2 


259.9 


259.9 


m eff (MeV) 


1072.4 


1084.8 


1063.4 


1073.1 


1068.6 


1070.6 


1069.8 


U (MeV) 


97.24 


88.88 


93.70 


91.44 


92.72 


92.13 


92.42 



Next, we investigate the A Am dependence of the order parameters. The ^i A ux depen- 
dence of < ipip > 1//3 and for iV=10 are summarized in Table ITU and Fig. |H1 The chiral 
condensate is quite stable to At. We show the A AVX dependence of the pion decay con- 
stants in two approximation schemes. In 'Scheme A' we neglect Y(p) and G(p) in ()2-38|) . 
'Scheme B' denotes the solution of (j2-48|) - (j2-52J) . The yl A ux dependence of n decay constant 
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is considerably decreased by the momentum scale expansion, 20% to 4%. 0(9°) quantity 
< ipip > 1//3 and m e ff are stable to the change of A AUX in this approximation. The stabilities 
of these quantities are maintained at finite temperature. For T « 120 MeV and N = Q, A Avyi 
dependence of < tpip > 1//3 was estimated 0.5% in the region At = [0.5, 1.0]. Therefore the 
critical temperature is also stable. Hereafter we choose At = 1.0 and N = 6. 

300 



250 
■^200 



-3- 

1-3-150 

V 



100 



50 



- i<!M>r 

• f„ (Scheme A) 
f, (Scheme B) 



0.8 0.9 



1.0 



0.4 0.5 0.6 0.7 

A 

Fig. 3. yl A ux dependence of the ir decay constant and the chiral condensate for the truncation 
N = 10. 



Table II. A AUX dependence of the order parameters. 



At 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


| < 4,4, > |V3 (MeV) 
/; [Scheme A] (MeV) 
/^[Scheme B](MeV) 


260.2 

130.47 

94.49 


260.1 

131.94 

96.65 


260.3 

130.17 

96.35 


260.3 

124.67 

95.08 


259.9 

117.59 

93.50 


260.0 

111.69 

93.21 


260.2 

105.30 

92.42 



3.2. Chiral phase transition at high temperature 

Let us analyze the chiral phase transition at high temperature. We first solve the ERG 
equation for the potential V(S) to estimate the critical temperature and the chiral conden- 
sate. The ERG equation for V(S) is already derived in ()2-37|) . Expanding the potential as: 
V = Vq + mS — g a S 2 /2 + GqS 3 /3 + ■ • ■ , and substituting to the ERG equation we have 

^o = C°ln(^ 2 + m 2 ), (3-5) 

A d , v 

A ~n m = -^r—zii ( 3 ' 6 ) 



alA A 2 + m 2 

4(-!) = ra(< 2 + 2m ^°-;ra>> (3 ' 7) 
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d (G 6 \ 2G 6 e , 2m f 3 x 8 m 3 C 3 1 

^dZ ItJ " - ^T^ + + » 2 ) 2 \ c +2mGeC ~3^T^J' (3 ' 8) 

^V-, (3-9) 



where denotes Matsubara frequency sums i.e. 

= ^NJffY! y/tf=d* (g c + ^trD- l /16N^N f ) 1 , (3-10) 



where iV/ = iV c = 3. 

The temperature dependence of the potential V are shown in Fig0] The order of phase 
transition is second order or very weak first order, and the critical temperature in our 
approximation is found to be 124.1 MeV. This result should be compared with T c = 169 
MeV given by Y. Taniguchi and Y. Yoshida, 10 - 1 T c = 166 MeV given by M. Harada and A. 
Shibata 11 -* and/or T c = 129 MeV proposed by O. Kiriyama, M. Maruyama and F. Takagi. 12 ^ 
Their analyses do not consider the gluon thermal mass. In Ref.13), the effect of the Debye 
screening at finite temperature and/or density is studied. The thermal mass is included in the 
static limit of the HTL approximation, thus the magnetic screening is ignored. T c gets 147 
MeV by this approximation. These works adopt the same running coupling constant g(p, k) 
at T = 0. Furthermore, the angular dependent part of the argument of g(p, k) = g((p + k) 2 ) 
was neglected. In our analysis, the higher order terms with respect to \p\ are not introduced 
at present. We consider the leading part of derivative expansion, and the argument of g 
becomes A in our analysis. At T = 0, this leading approximation scheme is identical to the 
prescription g(p, k) = g(max(j>, k)) in SDEs. 15 ^ 

Next, we estimate the pion decay constants at finite temperature. At finite temperature, 
there are two distinct 7r decay constants: the temporal one f n £, and the spatial one f n T- 
They are defined by 

j!-U L d°7T + --- } (3-11) 
ji~fnrd i TT + ---. (3-12) 

We can calculate these in terms of the wave function renormalization factor of meson field*- 1 

U L = y/ysZ^ao, U T = V^W^ixW (3-13) 

The ERG flow equations for Z^l and Z^t are already derived in ()2-48|) . Solving those 
equations, and substituting the solutions to (j3-13jl . we gain f n i and f n T- 
*) We renormalize the pion field as: 7Tr = Z^tt. 
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Fig. 4. Temperature dependence of the effective potential. The gluon thermal mass is incorpo- 
rated. 
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Fig. 5. Temperature dependence of the order parameters. 



The chiral condensate < ipip > can be evaluated by the common manner with the zero 
temperature case. The temperature dependence of the order parameters are shown in Fig. |3J 
| < tjjip > I 1 / 3 / f^{T = 0) keeps almost constant value 2.8 at low temperature T / f^iT = 0) < 
0.7 and decreases rapidly at Tj f w (T = 0) > 1.2. Finally, | < ipip > I 1 / 3 vanishes at T = T c . 
The transverse and longitudinal pion decay constants: f n x and f^L take almost same value 
at low temperature. The difference between them is less than 1% for T/f n (T = 0) < 0.4. 
However, for high temperature over T/f n (T = 0) > 0.8, the difference of two decay constants 
gradually becomes visible. Here, we can consider pion velocity in the thermal bath. The 
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pion momentum p and the pion energy po satisfy the dispersion relation: 

Vl = f%| a . (3-14) 

and pion's squared velocity v 2 is defined from the ratio of pion wave function renormalization 
factors, or decay constants: 

v 2 = Z^t/Z^l = f^r/ f vL . (3-15) 

The temperature dependence of the pion velocity v is drawn in FigEl It shows that pion 
travels at the velocity of light in the low temperature region, and that it slows down slightly 
as temperature rises. At T = T c , the pion velocity is slowed down to v/c « 0.95. 
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a. 
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Fig. 6. Temperature dependence of the pion velocity. 

Next, we want to know the critical exponents from the temperature dependence of the 
order parameters. From the construction in sec|21 our effective potential V is invariant 
under the transformation: ipip — > —ipip,mo — > — mo,o"o — > — ctq. Thus, the chiral condensate 
in ()3-3j) is considered to behave as: 

<W>=a f(<j 2 (T)). (3-16) 

< ipip > and o"o vanish at T = T c , however, /(0) keeps some finite value. Therefore, the 
critical exponent of the chiral condensate is derived by 

| < W > (T)| ~ Ca Q (T) ~ |T - T c \", (3-17) 

where C is some constant. The critical exponents of f n i and f^T are estimated in the same 
way Z nL and Z^t in (|3T3|) never go to zero when T — ► T c , so the exponents are determined 
by the temperature dependence of cr (T), same as in (|3T7|) . The critical exponents for 
the chiral condensate < ipip > and decay constants f n L and f n T should be universal. The 
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logarithms of | < ipip > j 1 ^ 3 //^ and |T — T c \/f n are plotted in FigEI We draw the line, 
log | < ipip > | = ulog \T C — T\ + C , where v and C are determined by the ordinary least 
squares, using the black points in Figd We can find the critical exponent: v = 0.48. 



1.0 - 
0.8 - 




Log U-Tcl/f* 

Fig. 7. The relation of log |T — T c \/f n and log | < ipip > l 1 ^ 3 //^- The black points are the datum 
used in the linear regression, and the white ones are not used. 

At zero temperature, it is verified that physical quantities do not depend too much on 
the choice of the IR regularization of the running coupling constant. 7 ^ So we similarly need 
to examine the Ai dependence of T c to guarantee the confidence of our results. In Ref.13), 
tf = ln(yli/ Aqcb) 2 dependence of the numerical results for the QCD phase structure is also 
discussed. It is shown that T c at \x = are stable for tf — 0.6, 0.5, and 0.4. In our calculation 
we find T c = 120 MeV, 124 MeV, and 129 MeV for t f = 0.6, 0.5, and 0.4 respectively 

3.3. Approximation scheme for the thermal mass 

So far, we carry out the calculations at finite temperature introducing the gluon thermal 
mass in the HTL approximation. This is the most striking difference from many previous 
works in the SDEs. There are several attempts to introduce the gluon thermal mass. We 
finally carry out the calculation in other approximations, and show how the difference is 
reflected to the result. In Table LTTTl we compare our result with other approximation schemes 
for the thermal screening effect. 

The values 'HTL approx' in Table LTTTl are our results and the case labeled 'po/|p| leading' 
is the well known form of the Landau damping, 24 ) 

m l = \ g ^{N c+ l -N f ^ rnl=\ml^. (3-18) 

This approximation has been applied to investigate the color superconductivity at high 
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Table III. The critical temperature and the critical exponent in several approximations. 



Scheme 


T c (MeV) 


V 


HTL approx. 


124.1 


0.48 


\Po\/\p\ leading 


98.6 


0.47 


Static limit of HTL 


154.3 


0.50 


Neglected 


171.9 


0.50 



density. 30 - ) Third one is the static limit of the HTL approximation, 8 )' 13 ) i.e. 

ml = \g 2 T 2 (n c + l -N^j , m 2 T = 0. (3-19) 

In the last case, the screening effect are dropped: 

ml = m 2 T = 0. (3-20) 

In Ref.ll) and Ref.12), the chiral phase transition at finite temperature/density has discussed 
in this approximation. T c in the second case is the lowest in our results. It is supposed that 
the factor |po|/|p| in ()3T8|) causes large uit for small \p\ (or large po) and strong screening 
effects. The critical exponents in the third case and the last case are found to be v = 0.50. It 
seems that the critical exponent, at least in our approximation, depends on the treatment of 
the magnetic mass. The deviation from our result are at most 4% for the critical exponent, 
while 20 ~ 40% for the critical temperature. The temperature dependences of the order 
parameters for each method are shown below. In FigJHl the temperature dependences of 
| < ijjip > I 1 / 3 is plotted. It tells us that the difference by approximation for the thermal 
mass appears at T/f n (T = 0) > 0.4. In each case | < ipip > | 1//3 begins to decrease at 
T/f^T = 0) « T c /f n (T = 0) - 0.4, and rapidly goes to zero over T/f^T = 0) « Tjf^T = 
0) -0.1. 

In FigEl we plot f n (T)/f n (T = 0), where f n (T) is so-called 'space-time averaged' decay 
constant, 11 ^ defined by /^(T) = (f 2 L + ^fnTf-KL) 1 ^ 2 /^- It tells that the effects of the thermal 
mass appear at Tjj^iT = 0) > 0.4. There is a definite peak in the 'Neglected' case. In 
SDEs, the temperature dependence of f^iT) is estimated, not that of f n L and f n T one by 
one. 11 -* In Ref.ll), the dependence is calculated by neglecting the thermal mass and there 
is a peak at T/f n (T = 0) « 1.2. 

The quantitative difference becomes evidently for the pion velocity v. In FiglTIH the 
difference of how v is slowed down is plotted, v/c is finally reduced to 0.92,0.92,0.95, and 
0.95 for 'Neglected', 'Static limit', 'HTL', and '|po|/|p| leading' respectively at each T c . v c 
at T = T c hardly depend on T c . It tells that v/c at T c depends obviously whether magnetic 
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mass is introduced or not. There are peaks in the 'Neglected' case and 'Static limit' case, 
and in the former case, v/c is increased to about 1.00 again near T/ f w (T = 0) = 0.9. 




0.0 0.5 1.0 1.5 2.0 

T/f„(T=0) 

Fig. 8. Temperature dependence of | < Y>?/> > | for each scheme. 



1.2 




T/f„(T=0) 

Fig. 9. Temperature dependence of fir(T)/f w for each scheme. 

§4. Summary and Discussions 

In this paper we analyzed the chiral phase structure of finite temperature QCD for 3- 
flavors by solving the sharp cutoff ERG flow equation. We derived the equation in the 
ladder approximation. The auxiliary fields a and 7r were introduced at the scale /Iauxj which 
is larger than /1 QC d- It was shown that the dependence of f n on T AUX was considerably 
improved by the momentum scale expansion. The gluon thermal mass is incorporated by 
the HTL approximation. As a result the critical temperature T c has become about 124 MeV, 
which is considerably lower value than the one calculated by neglecting the thermal mass. We 
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Fig. 10. Temperature dependence of pion velocity v for each scheme. 

also make certain that T c and the temperature dependence of the order parameters depend 
on the approximation scheme for the thermal mass. Our result shows that the screening 
effect is an important factor when the chiral phase structure of finite temperature QCD is 
discussed. 

In this paper, we replaced the gauge sector with the momentum dependent four fermi 
operator induced by one gluon exchange. Of course, the gauge dynamics in finite temperature 
QCD can be incorporated solely by the ERG. Then we have to employ the gauge invariant 
ERG 21 - 1 or MSTI. 20 ) As mentioned in ^21 both attempts will demand too much effort of 
human as well as machines. 

Including the instanton effect is the one of the interesting extensions of this work. There 
are various possible scenarios about the connection between U(1)a symmetry and the chiral 
phase transition. 31 - ) However, in the ERG framework U(l) symmetry breaking can not 
be incorporated straightforwardly because of the IR cutoff. 32 -* To study the effect of the 
anomaly, we have to add the instanton induced multi fermi term to at the initial scale. 

The phase transition seems to be second order, or very weak first order. We also studied 
the temperature dependence of the order parameters for each approximation scheme for 
the thermal mass. In the HTL, the chiral condensate and two 7r decay constants f n L and 
f-nT vanished for T « 124 MeV (T/f n (T = 0) ~ 1-34), that is the critical temperature. 
When T = 0, f W L and f n T have the same value. However, as temperature rises higher than 
T/f n (T = 0) ~ 0.8 and approaches T c , the difference becomes evident. The pion velocity v 
is determined from the ratio of the decay constants, and it decreases slightly as temperature 
approaches T c . 

It is one of the attractive subjects to apply our method to the case of finite density. 
In the passed few years, it has been clarified that the vacuum of hot/dense QCD has rich 



23 



structure, e.g. color super conducting phase, color-flavor locking phase, and so on. 33 -* The 
application of our method to the finite density case seems to be straightforward. 

Needless to say, HTL is based on perturbative QCD and it can not be said that this 
is truly proper method in the temperature region for the chiral phase transition, at most 
T « 100 MeV. However, this is just our first step of the exploration of the QCD phase 
structure by ERG. 
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Appendix A 

Thermal Threshold Functions 



In this appendix we give the explicit form of thermal threshold functions Iml,t) m (|2-48j) - 
(|2-52|) . First, we show what corresponds to the 0(p°) component and that of the spatial 
components. 



Jo = 2N f N c ^^^n (A2 ) Tv (A-l) 
-n 1 *-^ 1 [A z + m z ) 

n v ' 

n v ' 

n v ' 

The temporal threshold functions are a little complicated because of discretized expansion 
(gHH) and l[737j) . Iil and I 2L are given by 



hL = N f N c ^.^^^n (2 



(A 2 + m 2 ) 2 2ttT 



+N f N c ArVA 2 - tt 2 T 2 K u (A-4) 



(A-6) 



24 



where J\, J r 2, K\ and K 2 are defined as: 

_3 4 1 

1 ~ /l 2 + m 2 A 2 + m 2 + AixTuj n + 4vr 2 T 2 + /l 2 + m 2 + 87rTw n + 16vr 2 T 2 ' ^ ' ' 

32 = A 2 + m 2 ~~ A 2 + m 2 + A<uTuo n + 4tt 2 T 2 + yl 2 + m 2 + %irTu n + 16vr 2 T 2 ' 
2 1 

Kl = /l 2 + m 2 ~ A 2 + m 2 + 8tt 2 T 2 ' ( 

^2 = --r^ — ^ + T^ / » 9m9 - (A-10) 

A 2 + m 2 A 2 + m 2 + 8ir 2 T 2 

At T = 0, IkL and J^t reduce to an identical form. They are given by 

A 2 1 

2tt 2 yl 2 + m 2 v 7 

yl 2 yl 2 + 3m 2 

A _ 3^(yl 2 + m 2 ) 2 ' (A ' 12) 
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